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ABSTRACT 

We present a preliminary study of correlations in rapidity between pairs of identified charged 
pions, kaons and protons using the entire SLD data sample of 550,000 hadronic de- 
cays. Short range charge correlations are observed between all combinations of these hadron 
species, indicating local conservation of quantum numbers and charge ordering in the jet 
fragmentation process. The rapidity range of this effect is found to be independent of parti- 
cle momentum. A strong long-range K^-K^ correlation is observed at high-momentum and 
weaker long-range tt'^-tt' , tt'^-K^ and p-X^ and pp correlations are observed in light flavor 
events, providing new information on leading particle production in li, d and s jets. The long- 
range correlations observed in cc and hb events are markedly different and consistent with 
expectations based on known decay properties of the leading heavy hadrons. In addition, the 
SLC electron beam polarization is used to tag the quark hemisphere in each event, allowing 
the first study of rapidities signed such that positive rapidity is along the quark rather than 
antiquark direction. Distributions of ordered differences in signed rapidity between pairs of 
particles provide a direct probe of quantum number ordering along the quark-antiquark axis 
and other new insights into the fragmentation process. 
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1. Introduction 



Correlations between particles produced in hadronic jets can be used to probe details 
of the jet fragmentation process. In e^e~ annihilations into hadrons, the total charge, 
strangeness, baryon number, etc. of the final state particles in each event must be zero, 
and it is interesting to ask whether the conservation of such quantum numbers is local, 
or is longer-range in character. For example, in the case of strangeness, one can ask 
whether a strange particle and a corresponding antistrange particle tend to be produced 
"close" to each other within the event, whether strange and antistrange particles are 
associated with the initial s and s quarks, respectively, in ss events (or with the initial u 
and u in uu events, etc.), or whether strange and antistrange particles are distributed 
randomly throughout the event. Similar questions can be posed for other relevant 
quantum numbers. Previous studies of differences in rapidity between associated 
identified hadrons have shown that the conservation of charge, strangeness and baryon 
number is predominantly local. Most fragmentation models implicitly implement this 
feature, and the form and range of such short-range correlations provide useful tests 
of these models. Short-range correlations can also arise from the decays of heavier 
hadrons, for example the decay p° tt+tt^ will produce opposite-charge pion pairs 
with a characteristic degree of locality. 

Long-range correlations between particles of opposite charge and strangeness in 
opposite jets of an event have also been observed 0. These can be understood in terms 
of leading particle production whereby the higher-rapidity tracks in each jet tend to 
carry the quantum numbers of the initial quark or antiquark. In the case of e^e~ ss 
events, the s and s quarks may hadronize, for example, into a high momentum K~ 
and , respectively, and there need be no other strange particles in the event. In 
uu and dd events, however, the locality of quantum number conservation implies that 
a high-momentum strange-antistrange pair must be produced in each jet, which will 
dilute any long-range correlation. Nevertheless, improved measurements of long-range 
correlations may provide better understanding of leading particle production. 

The rapidity of a particle is typically defined with an arbitrary sign. If a sign could 
be given to each measured rapidity such that, for example, positive (negative) rapidity 
corresponds to the initial quark (antiquark) direction, then one might probe more 
deeply into both leading and nonleading particle production. One could measure, for 
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example, the extent to which a leading particle has higher rapidity than its associated 
antiparticle, and the extent to which low-momentum particles in jets remember the 
initial quark/antiquark direction. 

In this paper we present a study of correlations in rapidity between identified 
charged pions, kaons and protons based on about 550,000 hadronic Z° decays recorded 
by the SLD experiment at the SLAC Linear Collider. Clean samples of identified par- 
ticles were obtained using the Cherenkov Ring Imaging Detector. The hadronic event 
sample was divided into samples enriched in light-fiavor (e"'~e^ uu, dd or ss), cc 
and hh events in order to study the effects of the decays of the leading bottom and 
charmed hadrons. In section 4 we present a study of short-range correlations between 
all pair-combinations of these hadron species as a function of hadron momentum. We 
quantify the range of each correlation, and compare with the predictions of the JET- 
SET fragmentation model. In section 5 we search for long-range correlations between 
these species. In section 6 we exploit the SLC electron beam polarization to identify 
the quark (vs. antiquark) hemisphere in each event and sign the rapidities such that 
particles in the quark-tagged hemisphere have positive rapidity. The signed rapidity 
distributions themselves provide new information on leading particle production, and 
ordered rapidity differences between particle pairs allow new and unique probes of the 
fragmentation process. 

2. Apparatus and Hadronic Event Selection 

A general description of the SLD can be found elsewhere The trigger and initial 
selection criteria for hadronic decays are described in Ref. [Q. This analysis used 
charged tracks measured in the Central Drift Chamber (CDC) and Vertex Detector 
(VXD) 1^, and identified using the Cherenkov Ring Imaging Detector (CRID) []^. 
Momentum measurement is provided by a uniform axial magnetic field of 0.6T. The 
CDC and VXD give a momentum resolution of (Jp^/pi_ = 0.01 © 0.0026p^, where p± is 
the track momentum transverse to the beam axis in GeV/c. One quarter of the data 
were taken with the original vertex detector (VXD2), and the rest with the upgraded 
detector (VXDS). In the plane normal to the beamline the centroid of the micron- 
sized SLC IP was reconstructed from tracks in sets of approximately thirty sequential 
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hadronic decays to a precision of a/p ~ 7 fim for the VXD2 data and ^3 /im 
for the VXDS data. Including the uncertainty on the IP position, the resolution on 
the charged track impact parameter (6) projected in the plane perpendicular to the 
beamline is as =ll©70/(psin^/^ 6^) for VXD2 and as =8©29/(psin^/^ 6^) for 
VXDS, where 6 is the track polar angle with respect to the beamline. The GRID 
comprises two radiator systems that allow the identification of charged pions with 
high efficiency and purity in the momentum range 0.S-S5 GeV/c, charged kaons in 
the ranges 0.75-6 GeV/c and 9-S5 GeV/c, and protons in the ranges 0.75-6 GeV/c 
and 10-46 GeV/c §]. The event thrust axis p was calculated using energy clusters 
measured in the Liquid Argon Galorimeter ||10|| . 

A set of cuts was applied to the data to select well-measured tracks and events 
well contained within the detector acceptance. Gharged tracks were required to have a 
distance of closest approach transverse to the beam axis within 5 cm, and within 10 cm 
along the axis from the measured IP, as well as | cos^| < 0.80, and p± > 0.15 GeV/c. 
Events were required to have a minimum of seven such tracks, a thrust axis polar 
angle w.r.t. the beamline, 6t, within |cos^r| < 0.71, and a charged visible energy 
Eyis of at least 20 GeV, which was calculated from the selected tracks assigned the 
charged pion mass. The efficiency for selecting a well-contained —>■ qq{g) event was 
estimated to be above 96% independent of quark flavor. The VXD, GDG and GRID 
were required to be operational, resulting in a selected sample of roughly 285,000 events, 
with an estimated non-hadronic background contribution of 0.10 ± 0.05% dominated 
by Z^ — s> r+r~ events. 

Samples of events enriched in light and b primary flavors were selected based on 
charged track impact parameters 6 with respect to the IP in the plane transverse to the 
beam [O. For each event we define ngig as the number of tracks with impact parameter 
greater than three times its estimated error, 6 > 3as- Events with Usig = were 
assigned to the light-flavor sample and those with n^jg > 4 were assigned to the 6-flavor 
sample; the remaining events were classified as a c-flavor sample. The light-, c- and 
6-tagged samples comprised 166,000, 82,000 and S7,000 events, respectively; selection 
efficiencies and sample purities were estimated from our Monte Garlo simulation and 
are listed in table ^ 

Separate samples of hemispheres enriched in quark and antiquark jets were selected 
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Efficiency for — > 


Purity of 






UM, dd^ ss 


cc 


bb 


uu, dd, ss 


cc 


bb 


light-tag 


0.846 


0.338 


0.034 


0.881 


0.106 


0.013 


c-tag 


0.153 


0.617 


0.401 


0.320 


0.388 


0.292 


6-tag 


0.001 


0.045 


0.565 


0.005 


0.064 


0.931 



Table 1: Tagging efficiencies for simulated events in the three flavor categories to be 
tagged as light, c or b. The three rightmost columns indicate the composition of each 
simulated tagged sample assuming SM relative flavor production. 



by exploiting the large electroweak forward-backward production asymmetry wrt the 
beam direction. The event thrust axis was used to approximate the initial qq axis 
and was signed such that its ^-component was positive, > 0. Events in the central 
region of the detector, where the production asymmetry is small, were removed by the 
requirement \tz\ > 0.2, leaving 235,000 events. The quark-tagged hemisphere in events 
with left-(right-)handed electron beam was defined to comprise the set of tracks with 
positive (negative) momentum projection along the signed thrust axis. The remaining 
tracks in each event were defined to be in the antiquark-tagged hemisphere. The sign 
and magnitude of the electron beam polarization were measured for every event. For 
the selected event sample, the average magnitude of the polarization was 0.73. Using 
this value and assuming Standard Model couplings at tree-level, the purity of the 
quark-tagged sample is 0.73. 

For the purpose of estimating the efficiency and purity of the event flavor tagging 
and the particle identification, we made use of a detailed Monte Carlo (MC) simulation 



of the detector. The JETSET 7.4 ||T^ event generator was used, with parameter 
values tuned to hadronic e"'"e~ annihilation data |13[, combined with a simulation of B- 
hadron decays tuned to T(45') data and a simulation of the SLD based on GEANT 
3.21 [|l5l. Inclusive distributions of single-particle and event-topology observables in 
hadronic events were found to be well described by the simulation [ffl. 
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3. Identified Particle Selection 



The identification of cliarged tracks as pions, kaons or protons using tlie GRID is 
described in detail in [§. For tliis analysis we used a relatively loose identification 
algorithm, since the presence of misidentified hadrons at the 10% level has little effect 
on the measured correlations. Tracks with poor GRID information or that were likely 
to have scattered or interacted before exiting the GRID were removed by requiring 
each track to have at least 40 GDG hits, at least one of which was at a radius of at 
least 92 cm, to extrapolate through an active region of the appropriate GRID radiator 
and through a live GRID TPG, and, in the case of the gas radiator, to have fewer than 
four saturated hits within the volume in which the gas ring is expected. Approximately 
85% of the tracks within the GRID acceptance satisfied these cuts. 



Tracks identified in the calorimeters as electrons or muons were rejected, and for 



the remaining tracks log-likelihoods p|, [T^ were calculated for each of the three charged 



hadron hypotheses i = tc, K and p, and for each of the liquid and gas radiators. A 
track was tagged as a hadron of species i by the liquid (gas) system if the log-likelihood 
for hypothesis i exceeded both of the other log-likelihoods by at least 5 (3) units. In 
addition, for those tracks with good information from both the liquid and gas systems, 
the liquid and gas log-likelihoods were added together, and a track was tagged by the 
combined system if the log-likelihood for hypothesis i exceeded both of the others by 
at least 3 units. A track was identified as an i-hadron if it was tagged as type i by 
any of the liquid, gas or combined systems, and it was not tagged as any other type by 
any other system. The efficiencies for identifying accepted tracks are similar to those 
given in for charged pions there is roughly constant efficiency of about 80% in 
the momentum range 0.5-25 GeV/c; charged kaons and protons have similar efficiency 
except for a dip in the range 5-10 GeV/c. The simulation was found to provide a good 
description of the momentum distributions of the identified hadrons. 

For each identified track the rapidity y = 0.51n((i? + p\\)/{E — p\\)) was calculated 
using the measured momentum and its projection p\\ along the thrust axis, and the 
appropriate hadron mass. The distributions of rapidity are shown in fig. |^ for each 
identified hadron species, along with prediction of the simulation. Note that the overall 
sign of the thrust axis vector, and therefore the sign of the rapidity, is arbitrary. 
These distributions are not fiat in the central region, but show structure due to the 
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Figure 1: Rapidity distributions for the identified pions, kaons and protons. Also 
shown are the prediction of the Monte Carlo simulation. 



momentum dependence of the GRID identification efficiency (see The simulation 
provides a good qualitative description of the rapidity distributions; the dependence of 
the identification efficiency on momentum has only modest effects on the correlation 
studies, as discussed below. 

The absolute value of the difference between the rapidities of each pair of identified 
particles was calculated, and the distribution of this quantity is shown in fig. |^ for 
each of the six possible pairs of hadron species. In each case the distribution for those 
pairs with opposite charge is shown as the solid histogram and that for pairs with the 
same charge is shown as the dashed histogram. 
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Figure 2: Rapidity difference distributions for opposite-charge (histograms) and 
charge (dashed histograms) pairs of identified pions, Icaons and protons. 



8 



4. Short-Range Correlations 



For every type of hadron pair in fig. ^ there is an excess of opposite-charge pairs over 
same-charge pairs at small values of the absolute rapidity difference |Ay| = \yi — y2\- 
We expect more opposite-charge than same-charge track pairs, as well as more KK 
than KK /KK pairs and more pp than pp/pp pairs, due to conservation of electric 
charge, strangeness and baryon number, respectively. In the case of KK (pp) pairs the 
size of the excess is sensitive to the relative fraction of strange-antistrange (baryon- 
antibaryon) pairs that are both charged rather than neutral; if e.g. p-ii were the only 
type of baryon-antibaryon pair produced, there would be no such excess. The fact that 
the excess of opposite-charge KK and pp pairs peaks at low values of |Ay| indicates 
that conservation of strangeness and baryon number, respectively, is local, as has been 
observed previously 

In order to study these short-range correlations in more detail, we assumed that the 
tracks in each same-charge pair are unassociated, and subtracted their |A?/| distribu- 
tions from those of the respective opposite-charge pairs. These differences are shown 
for the low |Ay| region in fig. ^, and are seen to differ significantly from each other in 
form. Monte Carlo studies indicate that these differences are not due to acceptance, 
momentum dependence of the particle identification efficiency or to background from 
misidentified particles. Also shown in fig. ^ are the predictions of the simulation for 
this difference. In all cases, the simulation gives a reasonable description of the form 
of the difference, although there are small differences in the details of the form for vrvr, 
vrp, KK and vrp pairs. The predictions for the amplitudes have ~5% normalization 
uncertainties (not shown) due to uncertainties in the particle identification efficiencies. 
The predicted amplitudes are thus consistent with the data, except that for the i^p 
correlation, which is low by about 40%. 

Pairs of tracks from the decays of resonances (e.g. p'' tt+tt^) contribute to varying 
degrees for all pair combinations; each decay mode gives a |A?/| difference distribution 
with a characteristic form, and collectively they have a substantial influence on the 
overall form of the distributions. Since the resonance production rates are adjustable 
in the simulation, we suspect that it is able to provide an adequate description of the 
data. 

Pairs in which one or both particles are misidentified are an important source of 
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background, and lead to a convergence in the forms of the distributions for different 
pair combinations. The simulated differences for those pairs in which at least one track 
is misidentified are also shown in fig. ^ and are seen to contribute a sizeable fraction 
of the observed differences in the cases of unlike particles {tcK, vrp and Kp). There are 
normalization uncertainties on these predictions of ~20% due to uncertainties in the 
misidentification rates. Excellent particle identification is required to demonstrate, as 
seen in the figure, that there are excesses at short range for true pairs of these three 
combinations, indicating local conservation of electric charge between these different 
particle species and suggesting that there is charge ordering among hadrons of all 
species in the fragmentation process. The prediction of the simulation for the observed 
Kp correlation is dominated by misidentification; the prediction for the true correlation 
is therefore low by a factor of ~3. 

We have studied these short-range correlations in 6 bins of the momentum of the 
heavier (higher momentum) particle in the case of irK, vrp and Kp (vrvr, KK and pp) 
pairs. We find the quality of the simulated description of the data to be independent 
of the momentum bin. In order to quantify the range of each correlation, we fitted 
an ad hoc function, the sum of two Gaussians, to each difference in each bin over the 
range < |Ay| < 3 units. The centers of both Gaussians were fixed to zero, and 
the amplitude (width) of the wider Gaussian was fixed to 0.4 (2.2) times that of the 
narrower Gaussian, leaving the amplitude and width of the narrower Gaussian as free 
parameters. We used the width as a measure of the range of the correlation. This 
function provided a reasonable qualitative description of both the data and simulation 
for all pair combinations in all momentum bins. 

The fitted widths are shown in fig. ^ as a function of momentum for each pair 
combination. If the fragmentation process is scale invariant, then we expect the range 
to be independent of momentum, except for biases introduced by the bin edges and the 
particle identification. In the case of vrvr, the range is constant within ten percent except 
for the highest momentum bin. Significant momentum dependence is observed for 7rp, 
KK, Kp and pp pairs, however this dependence is reproduced by the simulation, with 
the possible exception of differences at low momentum for vrp and Kp pairs, and at high 
momentum for KK pairs. In the simulation, there is a ~10% momentum dependence 
due to binning and particle decays. The large slopes and point-to-point structure 
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Figure 3: Differences between the \Ay\ distributions for opposite-charge and same- 
charge pairs in the small |A?/| region. The errors are statistical only. Also shown are 
the predictions of the Monte Carlo simulation for all detected pairs (dashed lines) and 
those pairs in which one or both particles were misidentified (dotted lines). There 
are normalization uncertainties of ~5% on the total predictions and ~20% on the 
component due to misidentified particles. 
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are due to the momentum dependence of the particle identification efficiencies. Thus 
the hypothesis of momentum-independence of the range of the short-range rapidity 
correlation for a given pair combination is consistent with the data, within the context 
of the JETSET model. 

The qualitative results presented above were found to be independent of the primary 
ffavor of the event, as were the quality of the predictions of the simulation and the 
feature of scale independence. This is expected, as the only difference at short range 
should be due to the decay properties of the leading B and D hadrons; we do observe 
some small differences between the ffavor-tagged samples in the measured ranges of 
the correlations, which are reproduced by the simulation. 



5. Long- Range Correlations 

We next searched for long-range correlations between all pair combinations. In fig. H, 
a difference between opposite-charge and same-charge pairs at high \Ay\ is visible only 
in the case of KK, and even here the background from uncorrelated pairs is dominant. 
Since we expect long-range correlations from leading particle production to be relatively 
more important at high momentum, we have studied these differences for pairs in which 
both tracks have momentum p > 9 GeV/c. The corresponding \Ay\ distributions for 
each of the six pair combinations are shown in fig. ^. For these high-momentum pairs, 
there is a separation between pairs in the same jet (|A?/| < 2.5) and those in opposite 
jets (I Ay I > 4). At short range there is a large excess of opposite-charge pairs over 
same-charge pairs for all pair combinations, confirming that locality holds even at the 
highest momenta. At long range, there are clear correlations for tctt and KK pairs, as 
well as significant correlations for vrp, Kp and pp pairs. 

These long-range correlations are found to be strongly fiavor-dependent. Figure |^ 
shows the differences between the opposite-charge and same-charge pairs in the high 
I Ay I region for each of the ffavor-tagged samples. The bb events are seen to contribute 
very little to the difference for any pair combination, primarily because there are rela- 
tively few tracks with such high momentum in these events. Light ffavor and cc events 
contribute roughly equally to the observed correlations for tctt, vrp and pp pairs; light- 
ffavor events dominate the Kp correlation. In the case of ttK pairs, there is a strong 
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Figure 5: Observed \Ay\ distributions for opposite-charge (histograms) and same- 
charge (dashed histograms) pairs in which both tracks have p > 9 GeV/c. 



14 



Pair 




Data 






Simulation 




Type 


uds-tag 


c-tag 
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■nK 


171±73 


-245±45 


15±25 


43±34 


-169±20 


20± 9 


7rp 


103±53 


93±28 


-5±15 


145±25 


44±13 


8± 5 


KK 


552±33 


192±20 


37±12 


449±16 


190±10 


28± 4 


Kp 


101±33 


1±20 


-11±10 


76±15 


1± 9 


-6± 3 


PP 


39±17 


20±10 


2±10 


18± 7 


6± 4 


1± 1 



Table 2: Observed differences between the numbers of opposite-charge and same- 
charge pairs in which both tracks had p > 9 GeV/c and \Ay\ > 4, in events in each 
of the three flavor-tagged samples. Also shown are the predictions of the Monte Carlo 
simulation. 

anticorrelation in cc events along with a correlation in light-flavor events, both of which 
were invisible in the flavor-inclusive sample. 

The integrated differences for \Ay\ > 4 are given for each flavor-tagged sample in 
table ^ The predictions of the simulation are also given and are generally consistent 
with the data. The simulation also describes the rapidity dependence in fig. ^ within 
the experimental uncertainties. 

6. Signed Rapidities and Correlations 

We next tagged the quark (vs. antiquark) direction in each hadronic event using the 
electron beam polarization for that event, exploiting the large forward-backward quark 
production asymmetry in Z° decays. If the beam was left- (right-) handed then the 
thrust axis was signed such that cosO^ was positive (negative). Events with | cos6'j'| < 
0.15 were removed, as the production asymmetry is small in this region. The probability 
to tag the quark direction correctly in these events was 73%. 

The rapidity of a particle with respect to the signed thrust axis is naturally signed 
such that positive rapidity corresponds to the hemisphere in the tagged direction of the 
initial quark, and negative rapidity corresponds to the tagged antiquark hemisphere. 
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Figure 6: Observed differences between the |A|/| distributions for opposite-charge and 
same-charge pairs in which both tracks have p > 9 GeV/c, for the hght-(dots), c- (open 
circles) and h- (triangles) tagged samples. 
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Figure 7: Distributions (top) of the rapidity with respect to the signed thrust axis 
for positively (histogram) and negatively (dashed histograms) charged kaons. The 
difference (bottom) between these two distributions compared with the prediction of 
the Monte Carlo simulation. 

The signed rapidity distributions for identified and K~ are shown in fig. |^. There 
is a clear difference between the two, with more K~ than in the quark hemisphere, 
as expected due to leading K~ produced in s quark jets The difference between 

the two distributions is also shown in the figure and is compared with the prediction 
of the simulation, which is consistent with the data. 

For pairs of identified particles, one can define an ordered rapidity difference. For 
particle-antiparticle pairs, we define Ay^ = y+ ~ U- as the difference between the 
signed rapidities of the positively charged particle and the negatively charged particle. 
In fig. Hwe show the distribution of Ay^ for tt+tt", K^K~ and pp pairs. Asymmetries 
in these distributions are indications of ordering along the event axis, and the differences 
between the positive and negative sides of these distributions are also shown. The 
negative difference at high \Ay^ | for the K^K' pairs can be attributed to the fact 
that leading kaons are produced predominantly in ss events. Similar but smaller effects 
for expected for vr+vr" and pp pairs, and the positive difference in the latter at \Ay\ ^ 4 
may be attributed to this effect. For tt+tt^ pairs we observe a large positive difference at 
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high |A|/"''~| rather than the expected small negative difference, which is due entirely 
to cc events (see below). The predictions of the simulation are also shown and are 
consistent with the data at high \Ay^^\. 

The positive difference in the lowest \Ay^ | bins for the pp pairs indicates that the 
baryon in an associated baryon-antibaryon pair follows the quark direction more closely 
than the antibaryon. This could be due to leading baryon production and/or to baryon 
number ordering along the entire fragmentation chain. We find a significant effect in all 
six of our momentum bins, and that the bulk of the observed difference occurs at low 
momentum (see fig. We therefore conclude that both of these effects contribute; 
this is the first direct observation of baryon number ordering along the entire chain. 
Figure ^ also shows that the effect at long range noted above is confined to high- 
momentum pairs, as expected if due to leading baryon production. The prediction of 
the simulation is low by a factor of two at low |Ay+~| in both fig. ^ and in the low 
Pmax side of fig. |; it is consistent with the high Pmax data. 

In fig. |TU| we show the difference for K^K~ pairs in two bins of Pmax for the light- 
and c-tagged samples separately. For large Pmax one can see differences in amplitude 
and in |A?/| between the contributions from these two flavor-tagged samples. At low 
Pmax there is a positive difference of three standard deviations at low | Ay| for the light- 
flavor sample, a possible indication of strangeness ordering along the quark- ant iquark 
axis for associated K^K~ pairs, similar to the baryon number ordering observed for 
pp pairs. Such an effect is expected to be diluted at high momentum by an effect of the 
opposite sign due to triplets of high-momentum kaons produced in ss events. We do 
not observe such a signal in the heavy flavor events, possibly due to limited statistics, 
dilution due to heavy hadron decays that include a and a K~ , and/or dilution due 
to the negative difference from opposite hemisphere pairs, that populate a lower |A?/| 
region in heavy-flavor events than in light-flavor events. 

In fig. [Tl| we show the difference for tt+tt^ pairs in two bins of Pmax for the light- 
flavor and cc samples separately. For large Pmax we observe a negative difference at 
high |A|/| in the light-flavor sample, as expected from leading pion production. In the 
cc events, there are positive differences in both bins of Pmax at high |A?/+^| and also 
in the low Pmax bin at low |Ay+^|, all of which can be attributed to pairs involving a 
TT"*" from the D meson decay and a vr" from the D meson decay. The predictions of the 
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Figure 8: Distributions (left) of the difference between the signed rapidities of pos- 
itively and negatively charged identified hadrons of the same type. The differences 
(right) between the right and left sides of the distributions, compared with the predic- 
tion of the Monte Carlo simulation. 
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Differences between the positive and negative sides of the distributions of 
low (left) and high (right) momentum of the higher-momentum track. 
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Figure 10: Differences between the positive and negative sides of the distribution of 
yK+ — Vk- for pairs produced at low (left) and high (right) momentum, in the light- 
(dots) and c-tagged (open circles) samples. 
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Figure 11: Differences between the positive and negative sides of the distribution of 
— y^- at low (left) and high (right) momentum, in the light- (dots) and c-tagged 
(open circles) samples. 

simulation are consistent with these data. 

For unlike particles, both opposite-charge and same-charge pairs may be of interest. 
We define the ordered difference as the rapidity of the heavier particle minus that of 
the lighter particle multiplied by the charge of the heavier particle. That is, Ay^^ = 
yK+ - Vtt- or y^- - yK+ , Ay++ = yK+ -y^+, Ay— = y^- -yK~, etc. In fig. |l| we show 
the distributions of Ay^ for opposite-charge pairs of each of the three combinations 
of unlike particles, as well as the sum of the Ay^^ and Ay distributions for the 
corresponding same-charge pairs. A significant negative asymmetry is observed for 
ttK pairs of both opposite- and same-charge at all Ay, which may simply be due the 
combination of leading kaons and randomly selected pions. A similar effect in Kp 
pairs can be attributed to leading kaons combined with randomly chosen protons. In 
this case the asymmetry is negative for opposite-charge pairs and positive for same- 
charge pairs, as expected given our sign convention. In the case of vrp pairs there is a 
small positive asymmetry for both opposite- and same-charge pairs. The asymmetries 
predicted by the simulation are also shown; they are consistent with the pK~^ data, but 
overestimate the magnitudes of the effects for pK~, K^tt~ and K^t[^ pairs slightly. 
The latter two differences are expected in light of the absence of a t^K correlation in the 
simulation of light-flavor events noted above. A small negative asymmetry is predicted 
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for 7rp pairs, which is inconsistent with the observed positive asymmetry. 

7. Conclusion 

We have presented a prehminary study of rapidity differences between pairs of iden- 
tified charged pions, kaons and protons in hght-flavor hadronic decays. The SLD 
Cherenkov Ring Imaging Detector was used to select clean samples of identified charged 
hadrons, and the Vertex Detector was used to investigate the flavor dependence of the 
results. 

We observe excesses of opposite-charge over same-charge pairs for all pair com- 
binations at low values of the absolute rapidity difference, indicating that there is a 
high degree of local conservation of baryon number, strangeness and electric charge in 
the fragmentation process. The predictions of the JETSET fragmentation model are 
found to provide a qualitative description of the data, although they fail to describe the 
forms of some of the correlations in detail. The range of these short-range correlations 
was studied as a function of momentum; the variations found were reproduced by the 
simulation, verifying within the context of the JETSET model that the correlations 
are scale invariant. 

We observe a large excess of high momentum K'^K~ pairs over same-charge kaon 
pairs at large values of the absolute rapidity difference, and the effect is larger for higher 
momenta, as expected from leading kaon production in ss events. Weaker correlations 
are observed for protons, indicating that events with a leading baryon in one jet, a 
leading antibaryon in the other and no additional baryons do not contribute signifi- 
cantly to baryon production in e+e^ annihilations. Significant long-range correlations 
are observed between opposite-charge pairs of all combinations in light-flavor events; 
in cc events, long-range tttt, Trp and KK correlations are observed, along with a strong 
TiK anticorrelation. The simulation provides a good description of the data in general, 
but does not predict the long-range nK correlation in light-flavor events. 

We have studied distributions of rapidity signed so that positive rapidity corre- 
sponds to the quark (rather than antiquark) direction. Differences between signed ra- 
pidity distributions for positive and negative hadrons of all three species are observed, 
giving further evidence for leading production of charged pions kaons and protons. The 
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Figure 12: Distributions (left) of the difference between the signed rapidities of pos- 
itively and negatively charged identified hadrons of different types. The differences 
(right) between the right and left sides of the distributions, compared with the predic- 
tions of the Monte Carlo simulation. 
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distribution of the difference between tfie signed rapidities of and K" shows a large 
asymmetry at large values of the absolute rapidity difference, a direct indication that 
the long-range correlated KK pairs are dominated by ss events. A similar but smaller 
difference for tt+tt" pairs indicates roughly equal production of leading pions in uu and 
dd events. There is a large asymmetry at small rapidity difference for pp pairs, a clear 
indication of ordering of baryons along the event axis. A similar effect is observed for 
K'^K~ pairs at low momentum in light-flavor events. 
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